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Abstract 

The topological term in the effective action for the electrically neu- 
tral BCS system is discussed. It is applied for the calculation of the 
transverse force acting on the vortex in the limit of the smooth vor- 
tex core and vanishing interlevel distance in the vortex core. The 
controversy between the topological terms in cond-mat /9703124 and 
|cond-mat / 941 1047 is resolved. 
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1 Introduction 



There are many papers devoted to the effective action for the BCS system. 
They differ by the topological terms: while in most of papers the topological 
term is n<j>/2, where n is the particle density and is the phase of the 
order parameter (see eg |]]), in the several other papers it is (n — Co)<f>/2, 
where Cq = pp/2>ir 2 is the particle density in the normal metal (see eg H). 
From these different topological actions different predictions are made for 
the transverse force, acting on the moving vortex. Here we discuss this 
controversy. 

Actually both actions give the correct hydrodynamic equations for the 
BCS system in the uncharged limit. This is because the difference between 
the topological actions is the total time derivative if the parameter Co is 
treated as the dynamical invariant. This would mean that from both actions 
one should come to the same prediction for the vortex dynamics. However the 
latter depends on the proper treatment of the time derivative of the action. 
On the other hand, in the presence of the vortex, the phase is not defined 
globally, which makes means that the topological action in terms of <fi is 
badly defined. That is why we propose another form of the topological term 
effective action, which does not contain the phase of the order parameter 
explicitly, and which follows from the gradient expansion of the effective 
action. 

This action allows us to find the transverse force acting on the vortex in 
the special limit case, when the bosonic effective action provides the complete 
description of the vortex dynamics, ie in the case when the quantization of 
the levels in the vortex core can be neglected and the vortex core is smooth 
on the scale of the coherence length. In this limit our result for the transverse 
force coincides with that obtained in microscopic calculations ||. 



2 Effective BCS action at T = 0. 

The calculated topological term in BCS action can be expressed in terms of 
the order parameter (gap function) A in the following form: 

1 r ,s f 1 , f°° , 9n fdAOA* 0AdA*\ 
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Here r is an auxiliary coordinate, it is introduced when the effective action 
TrlnG -1 is presented as Tr drGdrG^ 1 . The dependence on r is chosen 
in such a way that at r = the coupling and the gap function are absent, 
A(r, t, t)| t=0 = 0, while r = 1 corresponds to the physical 3+1 dimensions, ie 
A(r, t, t)| t=1 = A(r, t). In Eq.(2.1) n is the particle density, which depends 
on A. At t = one has n{r = 0) = p|,/37r 2 = C , which is the particle 
number density in the absence of the gap at the same chemical potential. 

The quantity n — Cq is small in the weak coupling limit of the BCS model, 
where it is determined by the particle-hole asymmetry, n — Cq ~ nA 2 /Ep -C 
n. In a smooth crossover from the BCS superconductivity to the condensate 
of the Cooper pairs, discussed in many papers (see recent papers H, H, §), the 
parameter Cq decreases and can become zero at some value of the coupling 
parameter g: Co(g) = at g > g c . In this case at g = g c one has the quantum 
(Lifshitz) transition at T = 0. This zero temperature transition definitely 
happens if the quasiparticle spectrum has nodes, eg in the case of the <i-wave 
Cooper pairing or in the p-wave state with the symmetry of the A-phase 
of superfluid 3 He. In this case at g — g c the spectrum of quasiparticles is 
reconstructed - the nodes disappear (see Sections 6.2 and 9.4 in 0). This 
is similar to the Lifshitz transition in metals, where the topology of Fermi 
surface changes. In what follows we assume that we are on the BCS (weak 
coupling) side, ie below the Lifshitz point, g < g c . 

In the BCS action the relevant variable is the the gap function A, which 
means that all other variables, including n are the functions of |A| 2 . The 
variation of the action over the order parameter is the surface integral in the 
5-dimensional space (r, r, t) and thus is expressed in terms of the physical 
coordinates (r, t): 

In principle from the variation of this topological action one can restore 
the action in the physical space-time, but the well defined action contains 
explicitely the parameter C = n(|A| 2 = 0) = p 3 f /3tt: 

1 r rco 

Stop = j J d 3 r J ^ dt (n - n(\A\ 2 = 0)) 
If A(r, t) is nowhere zero (this requirement excludes the case when vor- 
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tices are present) one can introduce the phase <fi of the order parameter: 

A = |A|e^ . (2.3) 

This allows to express the variation of the topological action in Eq.(2.2a) in 
terms of the canonically conjugated variables, n and (f>: 

where 

5n = wk) m?) ■ (2 ' 5) 
The action in Eq.(2.2b) can be also used in this case, it gives 

1 r*. 



Sto P = --J d'r dt (n - C )-£ . (2.46) 

This is what was obtained in . 

The conventional contribution to the action is 

S = Jd 3 rJ2 o dt Qmn s v 2 + e(|A| 2 )) , (2.6) 

where n s is the superffuid density, m is the mass of the fermion, v s = ^;V0 
is the superffuid velocity and e is the energy density, which depends on |A| 2 . 
Varying S — So + S top over (f> and |A| 2 one obtains the conventional hydro- 
dynamic equations (we neglected the nonlinear terms in this procedure) 

— + V-(n s v s )=0 , (2.7) 

d<P 9 ^/g(|A| 2 ) 

It is important here, that the hydrodynamic equations are general and 
do not contain the parameter Cq. This is because the term containing the 
factor Co in Eq.(2.4b) is the full derivative. In this sense there is no difference 
between the topological terms discussed in |J and . The difference becomes 
important when the zeroes in A appear and the phase of the order parameter 
is not defined globally any more. This is just the case of vortices. To find 
the vortex dynamics we must return to the action in Eq.(2.2a) for the order 
parameter, which does not contain the phase (f) explicitly. 



4 



3 Effective action for the vortex dynamics in 
continuous limit. 



The effective action S = S + S top can be applied for the derivation of the 
dynamics of the vortex line only under certain conditions. We assume that 
the description in terms of the order parameter is complete, ie there is no 
other degrees of freedom. All the fermionic degrees of freedom are assumed 
to be inegrated off when the effective action was obtained. In this integra- 
tion the fermionic spectrum was considered classically, ie the spectrum was 
taken as a function of the commuting spatial coordinate r and momentum p: 
E = y / £ 2 (p) + |A(r)| 2 . This means that if one applies this effective action 
to vortices one neglects the quantization of the fermions in the vortex back- 
ground. This is justified only when the distance between the energy levels 
are small compared to the width of the levels, ujqt <C 1. This represents the 
necessary condition for the application of the effective action to the vortex 
motion. Another important condition is that the core size of the vortex is 
to be large compared to the coherence length in order to neglect the higher 
order gradient terms in the action. 

Under these conditions we can show that the BCS effective action leads 
to the following motion equation for the vortex: 

j[NTi{n — C )z x v L = irNhn s z x v s0 . (3.1) 

Here is the vortex velocity with respect to the heat bath (normal com- 
ponent or crystal lattice), which is here assumed to be at rest; v s0 is the 
superfluid velocity of the external superflow; N is the winding number of the 
vortex. This coincides with the result of the microscopic calculations for the 
electrically neutral case in the limit of the large core size and in the regime 

UJ T < 1 §. 

To get the Eq.(3.1) from the BCS action let us introduce the vortex 
coordinate (/!;). For simplicity we consider the rectilinear vortex along the 
axis z, so the order parameter A depends only on r — r^(t) where both vectors 
are 2-dimensional: 

A(r,i) = A(r-r L (t)) . (3.2) 
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The variation of the topological action in Eq.(2.2) becomes 

1 r l3 dn fdAdA* 8AdA*\ r°° , ( dx L . . dy L 
dS top = — I d r — - — ( — — — — ) I dt\ -^-by L - bx L 



2i J <9(|A| 2 ) \dx dy dy dx ) J-oo \ dt dt J 

(3.3) 

For the axisymmetric vortex with 

A{r,cp) = a(r)e tNtp , (3.4) 

where r and (p are cylindrical coordinates, and is the winding number, one 
obtains 

■ . dn fdAdA* dAdA*\ 
d 6 r 



d(\A\ 2 ) \ dx dy dy dx J 

= 2ttNL / dr— = 27iNL(n(oo) - n(0)) = 27iNL(n - C ) , (3.5) 
Jo dr 

where L is the length of the straight vortex line. We took into account that 
on the vortex axis (r = 0) the order parameter is zero, a(0) = 0, and therefore 
n(0) = C . 

The same result can be obtained for any type of the vortex with the 
winding number N, because the volume integral can be transformed to the 
surface integral far from the vortex core, which is determined only by N and 
n(|A(oo)| 2 ) - n(\A\ 2 = 0) = n - C . As a result 

5S top = 7rhNL(n - Co) jT dt {^-Sy L - Sx L ^f\ . (3.6) 

The variation of this topological action over 8yl gives the force acting on the 
vortex, which is represented by the term on the lhs of Eq.(3.1). 
Note that the expression for the variation of Eq.(2.2a): 

SS *>-L[*r ^L^VA'-VA^) (3.7) 



5r L 2i J d(\A\ 2 ) \ dt dt 

transforms to the Eq.(10) in if one neglects the coordinate dependence 
of dn/d(\A\ 2 ). The latter is justified only in the weak coupling limit, g 
1, explored in ||, while our approach is not limited by the weak coupling 
assumption: we assume only that g < g c . 
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In || the Eq.(lO) of |3j is generalized to the case of the anisotropic pairing 
with the momentum dependent gap function, A(p, r), but still in the weak 
coupling limit. We can generalize this to the arbitrary coupling strength 
introducing the momentum dependent particle distribution function 




d 3 p 
(2n) 



SS top = -i I — - I dtx 
a„(p,r) /-M(p,r) 8A-(p,r)\ _ 



n(p, r) = - | 1 - , = , Kv [ , = | , n(r) = 2 / j^n(p, r) . (3.8) 



Then 



<9(|A(p,r)| 2 ) V dt ' ^' ' dt 

Applying the same procedure as before one obtains the Eq.(3.6) for any 
(smooth) structure of the vortex core and for arbitrary anisotropic pairing 
state. I am indebted to N.B. Kopnin for this remark |5J. 

The rhs term in Eq.(3.1) is obtained from the kinetic energy term in So 
in Eq.(2.6). In the presence of the external superflow v s0 the relevant term 
is 

v s0 ■ / d 3 r n s v s = v s0 ■ P . (3.10) 



The linear momentum P related to the vortex coordinate is 

P = 7chNn s Lz x r L . (3.11) 

The variation of v s0 • P over gives the second term in Eq.(3.1). 

The equation (3.1) is not Galilean invariant. The Galilean invariance 
is restored by introducing the velocity of the normal component v n , which 
coincides with the velocity of crystal lattice in the case of superconductors, 
or with the heat bath of the normal excitations in the case of superffuids. 
This equation can be represented as the balance of three forces acting on the 
vortex [ H| : 



itTiNz x [n(v, - v L ) + C (v L - v n ) + (n - n s ){y n - v s )] = . (3.12) 
These are correspondingly the Magnus, the spectral flow and the Iordanskii 



forces in the terminology of Ref. [10|. 
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In the form of Eq.(3.1) the force balance was also obtained in ||1 1|| , again 
in the limit uoqt <C 1; the first discussion of the lhs term in Eq.(3.1) as orig- 
inating from the topological action was probably in [T2j (see the paragraph 
containing the Eq.(lO) in fT2||). 



4 Discussion. 

In conclusion, the effective topological actions discussed in j|] and f|| are 
equivalent to each other in the particular case when the phase of the conden- 
sate is defined globally, and thus in the absence of vortices. In this case the 
two topological actions differ by the total time derivative and thus lead to the 
same hydrodynamic equations for the conjugated variables, the particle den- 
sity n and the condensate phase <fi. However these actions cannot be applied 
when the zeroes in the order parameter are present. In the presence of zeroes 
the topological action in Eq.(2.2) (or in Eq.(3.9) for the more general case 
of anisotropic pairing) is to be used, which resolves between the discussed 
actions |lj and in favour of that in 

The Eq.(2.2) describes also the vortex dynamics in a special limit case. 
This dynamics agrees with the general phenomenological approach using the 
Poisson brackets scheme [jl3[ . In this approach the phenomenological param- 
eter Cq is introduced as the dynamical invariant of the system, which does 
not violate the general properties of the hydrodynamical Poisson brackets. 
The BCS effective action for the neutral BCS system in Eq.(2.2) gives the 
precise value of this phenomenological parameter, Cq = p^/3ir 2 , in the regime 
when the action is applicable for the vortex dynamics. The conditions of the 
applicability are ojqt <C 1 and R cor e 3> £o, where R cor e is the core radius. 
They correspond to the vortex with the smooth core. 

In the opposite limit lo$t = oo the effective action for the vortex dynamics 
in the BCS system can also be derived from the first principles exploring the 
Berry phase (see the recent paper [fLlfl). In this case the discrete spectrum of 



the electrons in the core is of main importance. The properties of the exact 
wave functions of the electron in the potential produced by the vortex lead to 
the zero value of the canonical momentum density at the origin, 5(0) = 0, in 



TM. This gives the zero value of the parameter Cq in Eq.(3.12) in complete 



agreement with the result of [10], obtained in the limit Uqt = oo. The finite 



life time r destroys the coherence of the wave functions and finally restores 
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the maximal value of Co in the limit ujqt <C 1. 

I would like to thank Jacek Dziarmaga, Nikolai Kopnin and Anne van 
Otterlo for stimulating discussions. 
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